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Abstract: In Atomic force microscope (AFM) examination of a vibrating microcantilever, 
the nonlinear tip-sample interaction would greatly influence the dynamics of the cantilever. 
In this paper, the nonlinear dynamics and chaos of a tip-sample dynamic system being run 
in the tapping mode (TM) were investigated by considering the effects of hydrodynamic 
loading and squeeze film damping. The microcantilever was modeled as a spring-mass- 
damping system and the interaction between the tip and the sample was described by the 
Lennard-Jones (LJ) potential. The fundamental frequency and quality factor were 
calculated from the transient oscillations of the microcantilever vibrating in air. Numerical 
simulations were carried out to study the coupled nonlinear dynamic system using the 
bifurcation diagram, Poincare maps, largest Lyapunov exponent, phase portraits and time 
histories. Results indicated the occurrence of periodic and chaotic motions and provided a 
comprehensive understanding of the hydrodynamic loading of microcantilevers. It was 
demonstrated that the coupled dynamic system will experience complex nonlinear 
oscillation as the system parameters change and the effect of squeeze film damping is not 
negligible on the micro-scale. 

Keywords: TM-AFM; microcantilever; squeeze film damping; Lennard-Jones (LJ) 
potential 
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1. Introduction 

Atomic force microscopy (AFM) has been developed to a nearly ubiquitous tool for studying 
physics, chemistry, biology, medicine and engineering at the nano-scale [1-5]. AFM could 
significantly impact many fabrication and manufacturing processes due to its advantages such as 3D 
topography of nano-fabrication and metrology for MEMS [2]. As a typical dynamic mode, the tapping 
mode (TM) is widely used in the operation of AFM where the cantilever is driven at a fixed frequency 
close or equal to the fundamental resonance frequency of vertical bending [4,6]; a schematic of the 
TM-AFM setup is shown in Figure 1. Moreover, the vibration amplitude of the cantilever is much 
bigger than the equilibrium separation between the tip and the sample. The TM-AFM has attracted 
extensive attention due to its ability to deal with compliant materials as well as to overcome adhesion 
forces. 

Figure 1. Schematic of the typical TM-AFM setup. 



Photodetector 




The inherently and highly nonlinear tip-sample interaction will give rise to complex dynamics of 
the cantilever in TM-AFM [7]. Nonlinearity is essential in understanding the dynamics of cantilevers 
as there are many nonlinear forces in TM-AFM, such as the attractive van der Waals forces, the short- 
range repulsive interactions, contact nonlinearities and capillary forces et al [8-10]. Ashhab et al [11] 
concluded that chaos in AFM depend on the damping, excitation and tip-sample distance, and 
suggested that a state feedback control can be used to eliminate the possibility of chaotic behavior. 
Using nonlinear analysis methods and numerical simulations. Basso et al [12] found that the chaotic 
behavior may occur via a cascade of period doubling bifurcations. In their studies [11,12], Melnikov 
theory was used to predict the existence of chaos in AFM. The nonlinear dynamics to frequency 
sweeps in TM-AFM was simulated using the van der Waals forces and the Derjaguin-Muller-Toporov 
(DMT) and Johnson-Kendall-Roberts (JKR) contact models [6,13]. Lee et al [6] carried out numerical 
analysis using modern continuation tools for computational nonlinear dynamics and bifurcation 
problems where the tip-surface interaction was represented by the van der Waals and DMT contact 
forces. Nonlinear hysteresis and jumps in the dynamic response were examined as the tip approaches 
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to or retracts from the sample at a fixed excitation frequency [14]. Zitzler et al [15] considered the 
infiuence of hysteretic capillary forces in TM-AFM and studied the effect of the relative humidity on 
the amplitude and phase of the cantilever oscillation. By using the forward-time simulation and 
numerical continuation techniques [15], Hashemi et al [10] investigated the nonlinear dynamics of a 
TM-AFM with tip-surface interactions, which include attractive, repulsive, and capillary forces. 
Hersam [3] conducted experiments to verify the importance of nonlinear dynamics in TM-AFM 
measurements. Rutzel et al [9] used the Lennard-Jones (LJ) potential to model the tip-surface 
interactions and carried out a comprehensive investigation to the nonlinear dynamics and stability of 
the TM-AFM, and the results showed that considering the LJ interaction potential in modeling the 
dynamics of AFM could improve the qualitative prediction of the real system response. 

Nonlinear dynamics and chaos of the cantilever for the TM-AFM remain a challenging and critical 
issue. After experimental observations. Couturier et al [16] mentioned that the motion of the 
cantilever could become chaotic under instability conditions. When the cantilever is driven close to the 
surface of the sample, the squeeze film between the cantilever and the sample surface contributes 
significantly to the damping and gives rise to the complicated nonlinear behavior [17,18]. When 
studying the frequency response of AFM cantilevers in liquid media contained in a commercial fiuid 
cell, Motamedi and Wood- Adams [19] found that such systems could exhibit complicated dynamics. 
However, a rational connection of the tip-sample-interaction and the nonlinear dynamics analysis of 
the cantilever where the coupled effects of squeeze film damping and hydrodynamic loading are 
considered has not been presented and addressed satisfactorily. To order to make the TM-AFM 
achieve good performance, it is necessary to identify and so as to eliminate the possible chaotic motion 
of the cantilever of the AFM. In this paper, the cantilever is modeled as a single spring-mass-damper 
system and a nonlinear dynamic model is developed to study the cantilever-sample interaction by 
using the LJ potential including the long-range attractive forces and short-range repulsive forces. A 
comprehensive investigation of the nonlinear dynamics and chaos of the TM-AFM is carried out. 

The rest of the paper is organized as follows. Section 2 describes the mathematic model of the 
cantilever vibrating in air considering the effects of hydrodynamic loading and squeeze film damping 
and the physical model of the cantilever-sample interaction is established in Section 3 using the LJ 
potential. Numerical results and discussions of the quality factor and resonant frequency of the 
frequency response, and nonlinear chaos and bifurcation of the dynamic TM-AFM are presented in 
Section 4. Finally, we end the paper with our conclusions in Section 5 . 

2. Micro- Cantilever Vibrating in Air 

2.1. Hydrodynamic Loading Effect 

The gas flow around the cantilever is assumed to be incompressible and the Navier-Stokes equation 
is given by [20]: 

where p , y , p and // are the density, velocity, pressure and viscosity of the gas, respectively. 
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It has been found that the wavelength of vibration greatly exceeds the dominant length scale in the 
flow [20,21]. The nominal width B is the dominant length scale in the gas flow and the appropriate 
RejTiolds number Rg can be expressed as [21]: 

_ pcoB^ 

where is a characteristic radial vibration frequency and Re is a normaUzed Reynolds number which 
indicates the importance of viscous forces relative to inertial forces in the gas fluid. It can be found 
that the gas could be considered to be non- viscous in the limit as Re oo, and the non-viscous gas 
model is applicable for practical cases where Re » 1 [21]. However, if co is close to the resonant 
frequency of the cantilever in vacuum, a reduction in dimensions of the beam will result in a reduction 
in Re. 

The surrounding viscous fluid medium plays an important role on the dynamics of the cantilever in 
AFM, Sader et al [21,22] studied the effect of viscous fluid medium on the AFM cantilever and found 
that the shift in resonant frequency of the cantilever from vacuum to fluid (gas or liquid) was strongly 
dependent on both the density and viscosity of the fluid. For a rectangular cantilever beam [21], when 
the quality factor of the fundamental mode of the cantilever in gas exceeds 1, which is typically 
satisfied when the cantilever is placed in air, the relationship between the vacuum resonant frequency 
coyac and the resonant frequency in gas cogas can be written as [20]: 



CO 

gas 



npB 



V 



(3) 



where pc and H are the density and thickness of the cantilever, respectively, and the natural scaling 
parameter is n = -^, which is defined as the ratio of the added mass of the gas to the mass of 

the cantilever. 

To study the effect of hydrodynamic loading on the dynamics of the cantilever, the hydrodynamic 
functions r(a;), which represents the real and imaginary pressure of the surrounding on the cantilever 
in two dimensions, is given by [20,21]: 

r = r,+7T, (4) 

where r, =a, and r, =-^+^, and ai = 1.0553, a2 = 3.7997, bi = 3.8018, and \>2 = 2.7364 are 



selected from [20]. 

The quality factor for the bending modes can be written as [20]: 

T.(co ) r.(co ) 

i ^ gas ' I \ gas / 

For small amplitude of the normal vibration of the cantilever, the local interaction stiffness in the 
vertical direction kn and damping coefficient C„ is given by [23]: 



^C0S^3-1 

J 



(6) 
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And: 



(7) 



where A^o is the free amplitude in the normal direction, is the measured amplitude, cp is the 
measured phase of the cantilever, co is the drive frequency, and ki is the normal bending stiffness and 



EBH' 



, in which E, B, H and L are the elastic modulus, width, thickness and length of the 



cantilever, respectively. When the cantilever in AFM operates far below the resonance, the phase angle 
is close to zero, then kn and C„ become: 



k^^kAK-\) C =0 

n L \ / r> n 



(8) 



where k = A^^l A^, 



2.2. Squeeze Film Damping Effect 



When the cantilever is driven close to the sample, the squeeze film between the cantilever and the 
surface of the sample causes significant damping except for the fluid dissipation at the edges and 
above the cantilever. For the damping problems encountered at micro-scale [24], the squeeze film 
damping can be expressed by the nonlinear Reynolds equation: 

d{hp) 



_d_ 



M 5^ 



+ - 



¥ dp 

P — 

M 5^, 



= 12- 



dt 



(9) 



Under the basic assumptions including the gas in the gap has been regarded as a continuum and the 
gas undergoes an isothermal process, the squeeze film damping can be modeled as: 



Pah' 



■eff 



Pa 



+ - 



Pa 



dt 



+ - 



dt 



Pa 



(10) 



where the pressure p is small compared to the ambient pressure pa (p^ =1.013xlO^Pa ), is the 
effective gas viscosity, and h is the gap between the cantilever and the surface of the sample. At micro- 
scale, when the ratio between the mean free path of the air X and the film thickness h, i.e., Knudsen 
number Kn = X/h, is not small, the no-slip boundary condition at the interface may be inadequate. 
Considering the effect of slip fiow, the effective viscosity jUe/f can be expressed as [25]: 



^eff 



1 + 9.638^, 



1.159 



(11) 



where ju is the viscosity coefficient at 1 atm. 

For the normal motion of parallel plates, h and jUef/avQ not functions of the position. Equation (10) 
can be simplified as: 



12//, 



P 



dp 



^ d 

H 



dp 



d(ph) 

dt 



(12) 
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Then, the damping pressure can be obtained by direct integration with the boundary conditions, i.e.: 



n 



dh 
dt 



(13) 



The damping force on the cantilever is: 

F^=f_',lp{^,t)Ld^ = 
and the coefficient of damping force is given by: 

C = 



dt 



^effB'L 



(14) 



(15) 



3. The Physical Model 

The AFM is composed of an elastic cantilever and the achievable sensitivity and resolution of AFM 
depend largely on the geometry of the cantilever [11,26]. Considering only the first vibration mode, 
the cantilever can be modeled as a simplified spring-mass-damping system, as shown in Figure 2. The 
tip is modeled as a sphere of radius i?, and the cantilever-sample distance is characterized by zq, which 
is the distance between the equilibrium position of the cantilever and the sample when only the gravity 
acts on it. The cantilever position is given by x measured from the equilibrium position. 

Figure 2. Schematic of the lumped spring-mass-damping model for the TM-AFM 
cantilever vibrating near a sample surface. 
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The LJ potential models the dispersive van der Waals forces as well as the short-range repulsive 
exchange interactions between two molecules [9]. The interaction between a cantilever tip and sample 
surface can be modeled as the interaction between a sphere and a flat surface. The tip-sample 
interaction is modeled by the LJ potential given by [9,1 1]: 

1260(zo+x) 6(zo+x) ^'^^ 
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where Ai and A2 are the Hamaker constants for the attractive and repulsive potentials, respectively. 
The Hamaker constants are defined as ^1 = Tip\p2C\ and A2 = npipici, in which pi and p2 are the 
densities of the two interaction components, and ci and C2 are the interaction constants, respectively. 

When the cantilever is driven close to the sample, there are three different kinds of forces including 
the spring force, the van deer Waals attractive force which is proportional to the inverse square power 
of the distance between the cantilever tip and the sample, and the repulsive force which is proportional 
to the inverse eighth power of the distance between the tip and the sample. The LJ force can be defined 
as the sum of the attractive and repulsive forces and is expressed as [9]: 

^ dU,j{x,z,) ^ A,R A,R 

e(x+z,) m{z,^xf 6{z,^xf ^^^^ 

During the AFM operates in the TM, a low-dimensional model reduction can provide an accurate 
description of the cantilever dynamics. As shown in Figure 2, the cantilever is driven by the harmonic 
driving force, the tip-sample interaction force Fu (LJ force) and the force due to squeeze film damping 
Fs, The governing equation of the motion of the cantilever can be determined by: 

mx + cx + kx + k^x^ =F^^(x,Zq) + F (x,x,Zo) + f^coscot (18) 

where x is the instantaneous displacement of the cantilever tip measured from the equilibrium tip 
position in the absence of external forces with positive values toward the sample surface, and x and 
X are the instantaneous velocity and acceleration of the cantilever tip,, m, k and c are the equivalent 
mass, spring stiffness and damping coefficients of the cantilever in air, respectively, and k ^k^+k^ 

and c = mco^ I Q + C^, in which co^ = -sjkjm is the first-order mode frequency, and kc is the nonlinear 

cubic stiffness. /o and co are the amplitude and angular frequency of the harmonic driving force. 

The damping force F^(x,x,Zq) due to the squeeze film on the cantilever can be obtained from 

equation (14) and is given by: 

F (x,x,Zo ) = c^x = -^x (19) 

(x + zj 

In order to facilitate the investigation of the qualitative behaviors of the dynamic system, the 
equilibrium distance variable Zs is defined and = (3/2)(2£))^^^ [1 1], in which D = A^R/(6k) . 

Introducing dimensionless variables: 

^ . X 4 z^ k 0 



co,Z^ 27 Z'k mco, 



A 



\A J 



(20) 



Then, the LJ force Flj and the squeeze film damping force Fs can be rewritten as: 

d 

Fu-~~ ^"""TIT ^ (21) 
{a + y) 30(a + y) 



And: 
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Fs-X- (22) 



Therefore, the dynamic equation of the system can be given by: 



y + Cy + y + Py 



r \ 



FcosDr-// ^- — -y 

(« + j) 



3 

J 



(23) 



{a + y) 30(a + 3;) 
where, F = Fq/s, k^yl s , ri = xl ^ ^ ^i^d e is a small perturbation. 

4. Results and Discussion 

This section aims at numerically investigating the characteristics and nonlinear dynamics of a TM- 
AFM cantilever-sample system driven by the harmonic excitation. The general properties of the 
cantilever and interaction properties with the respective sample are referred to [9], are listed in Table 
1. The 4th order Runge-Kutta method is used to integrate the set of Equation (23). A small integration 
step (27i/200) has to be chosen to ensure a stable solution and to avoid the numerical divergence at the 
points where derivatives of Fu and Fs are discontinuous. The effects of system parameters on the 
dynamic behavior of the cantilever vibrating system are investigated by using the bifurcation diagram, 
Poincare maps, largest Lyapunov exponent, phase portraits, time histories and amplitude spectrum. 

Table 1. Parameters of the silicon tip tapping and silicon sample used in the 
numerical simulations. 



Description 


Value 


Length 


449 |im 


Width 


46 |im 


Thickness 


1.7 |im 


Tip radius 


150 nm 


Material density 


2,330 kg/m^ 


Young's modulus 


176 GPa 


Bending stiffness 


0.11 Nm-^ 


First-order resonant frequency 


11.804 kHz 


Quality factor 


100 


Hamaker constant (Repulsive) 


1.3596 X 10-^^ J-m^ 


Hamaker constant (Attractive) 


1.865 X 10"^^ J 



4. 1. Effect of External Forcing Term F 



The external forcing term is one of the most important parameters affecting the dynamic 
characteristics of the TM-AFM tip-sample system. Figure 3 and Figure 4 show the bifurcation diagram 
and largest Lyapunov exponent map of the dynamic system where the amplitude of the external 
excitation is the control parameter and the small perturbation ^ = 0. 1 is added. 
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Figure 3. Bifurcation diagram of the amplitude of the external forcing term F. 
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Figure 4. Largest Lyapunov exponent map of the amplitude of the external forcing term Y. 
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It can be seen from Figure 3 that the system responses contain periodic and chaotic motions 
alternately at the interval of 0 < F< 80. When F = 6.8, the vibration amplitude of the cantilever is 
small, the motion is synchronous with period-eight (P-8), and eight points are correspondingly 
displayed in the Poincare map, as shown in Figure 5(a). With the increase of the amplitude of the 
forcing term F , the motion becomes synchronous with period-two (P-2) at F = 33.0, as illustrated in 
Figure 5(b). Moreover, the period- 1 motion with only one isolated point in Poincare map and one 
circle in phase portrait can be observed at F = 65.0 in Figure 5(c), and the corresponding largest 
Lyapunov exponent becomes negative, according to Figure 4. At F = 75.2, the chaos is shown in 
Figure 5(d), the strange attractor has a fractal structure in Poincare map, and it can be found in Figure 
4 that the corresponding largest Lyapunov exponent is positive. Therefore, as the amplitude of the 
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forcing term increases, the changes of the system responses are very complex, with alternative periodic 
and chaotic motions. 

Figure 5. The Poincare maps, phase portraits, amplitude spectrums and time histories of 
different F. 



(a) r = 6.8 



(b) r = 33.o 



(c) r = 65.0 



(d) r = 75.2 
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4.2. Effect of Squeeze Film Damping rj 



At the micro-scale, the squeeze film damping coefficient and ratio are the key parameters for the 
dynamic responses of the micro-devices. The larger the squeeze film damping is, the higher the noise 
level results. The ratio of the fundamental resonant frequency in gas to that in vacuum cogajcovac is 
numerically calculated from Equation (3). Figure 6 gives the ratio of resonant frequencies in vacuum 
and gas cOgaJcOvac as a function of the Reynolds number Re at different natural scaling parameter n . It 
could be found that cOgaJco^ac is increasing with the increase of Reynolds number Re and cOgaJcOvac 
decreases with the increase of the scaling parameter FL. Figure 7 gives the quality factor Qgas as a 
function of both natural scaling parameters. Re and FL, which can be obtained directly from Equation 
(5). It is indicated that the quality factor Qgas increases with the increase of the Reynolds number Re 
and the decrease of the natural scaling parameter FL. However, when the natural scaling parameter 11 
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tends to be a larger value, i. e., IT = 10, the quality factor Qgas has small change. For example, when IT 
changes from 10 to 100, the quality factor Qgas has very small change, as shown in Figure 7. 

Figure 6. The relationship between the vacuum resonant frequency co^ac and the resonant 
frequency in gas cogas (cogajcovac) as a function of the Reynolds number Re at different 
natural scaling parameter H. 
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Figure 7. The quality factor Qgas as a function of the Reynolds number Re for the 
fundamental mode at different natural scaling parameter H. 
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Microstructures undergoing motion transverse to a fixed plate exhibit damping effects that should 
be considered in the dynamics simulation. The heat transfer analogy is applied to obtain the damping 
and stiffness coefficients of the microstructures with PLANE 55 thermal elements in ANSYS 11.0 
[28]. The effective damping and stiffness coefficients are determined by the finite element thermal 
analogy approach and theoretic analyses at different operation frequencies and the results are listed in 
Table 2. It can be seen that the damping coefficient decreases tardily with the increase of the operation 
frequency and the stiffness coefficient increases fleetly at the same time with slip and without slip. 
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Meanwhile, the damping and stiffness coefficients with slip effect are smaller than those without slip. 
Figure 8 shows the pressure distribution of the damping component with slip at the resonant frequency. 
The pressure distribution is approximately parabola in the directions of length and width, and the peak 
appears at the center of the film. 



Table 2. Effective damping and stiffness coefficients of the squeeze film at different 
frequencies. 



Frequency 



Damping coefficient 



Stiffness coefficient 



PLANE55 



Analytic(slip 



PLANE55 



Analytic(slip 





No slip 


Slip 




)[27] 




No slip 


Slip 


)[27] 


1 


1.5199e-4 


1.3499e- 


4 


1.2529e 


-4 


9.2051e- 10 


7.2610e- 10 


6.2304e- 10 


1 000 


1.5199e-4 


1.3499e- 


4 


1.2529e 


-4 


9.205 le- 4 


7.2610e-4 


6.2304e - 4 


11 804 


1.5197e-4 


1.3497e- 


4 


1.2528e 


-4 


0.1282 


0.1012 


0.0868 


50 000 


1.5163e-4 


1.3473e- 


4 


1.2509e 


-4 


2.296 


1.8118 


1.555 


100 000 


1.5056e-4 


1.3398e- 


4 


1.2449e 


-4 


9.117 


7.206 


6.190 



Figure 8. Pressure distributions of the cantilever at the resonant frequency. 
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Figure 9. Bifurcation diagram of the squeeze film damping ratio 77 
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Figure 9 displays the bifurcation diagrams of squeeze film damping ratio rj in the range of 
0.05 < // < 0.1 for the coupling nonlinear dynamic system, from which it can be that the system 
response changes between periodic and chaotic motions alternately. Figure 10 shows the Poincare 
maps and phase plane portraits of different squeeze film damping ratio rj on the responses of the 
coupling system. The system response starts synchronous motion with period-4 at // = 0.0525 (shown 
in Figure 10a), and then becomes synchronous motion with period- 1 at rj = 0.0635(shown in Figure 
10b), and then leaves synchronous motion with period- 1 and enters chaotic motion at // = 0.07, which 
can be seen in Figure 10c. The strange attractor has a fractal structure and the corresponding largest 
Lyapunov exponent is positive. As indicated in Figure lOd, with the increase of squeeze film damping 
ratio, the system response becomes synchronous motion with period- 1 from chaotic motion. Therefore, 
the effect of squeeze film damping on the system response cannot be neglected for structures at the 
micro-scale. 

Figure 10. The Poincare maps and phase portraits of different squeeze film damping 
ratios //. 



(a) r] = 0.0525 



(b) T] = 0.0635 



(c) rj = 0.07 



(d) r] = 0.0835 
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4.3. Effect of material property parameter Z 



At micro-scale, the material properties of the AFM tip and sample play an important role on the 
surface force between them and, as a result, the dynamic response of the tip-sample system displays 
very rich nonlinear characteristics. Figure 11 is the bifurcation diagram of the material property 
parameter E for the coupling nonlinear dynamic system with different cubic stiffness ratios and 
squeeze film damping ratios. The system parameters are taken as follows: equilibrium parameter 
a= 1.2, excitation frequency ratio Q = 1, integration step numbers = 200 and the bifurcation step 
AS = 0.005 . 

Figure 11a shows the bifurcation diagram of the material property parameter in the range of 
0.3 < i7 < 0.42 with the cubic stiffness ratio y? = 0.3. The response of the coupled nonlinear system 
undergoes a complete process from chaotic motion through period-2, period- 1, period-4 and period-8 
motions to steady-state motion with period- 1 by the forms of period-doubling and anti period-doubling 
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bifurcations. At the interval of 0.42 < Z < 0.47, the system response alters between chaotic motion 
with long time and periodic motion with short time. When Z increases to i7 > 0.47, the system 
response becomes synchronous motion with period- 1. With the change of the cubic stiffness ratio of 
the cantilever tip from y? = 0.3 to y? = 0.5, the response of the coupled system undergoes the process 
of chaotic and periodic motions alternatively. It comes into period-2 motion from chaotic motion with 
anti period-doubling bifurcation, and then enters period- 1 motion in a large range of material property 
parameter > 0.39), as illustrated in Figure lib. In addition, it is found that the chaotic motion 
disappears at higher i7 with the increase of cubic stiffness ratio (from y? = 0.3 to y? = 0.5). Therefore, 
with the changes of the measured samples in experimental tests, the values of i7 vary accordingly, and 
as a result the response of the coupled system displays various nonlinear dynamic behaviors. 

Figure 11. Bifurcation diagram of th 1 cm e material property parameter i7 at different 
cubic stiffness ratios: (a) y? = 0.3; (b) y? = 0.5. 




Z Z 



Figure 12 is the bifurcation diagram of the cubic stiffness ratio y? on the response of TM-AFM tip- 
sample system at the interval of 0.3 < yff < 0.6 for various combinations of squeeze film damping 
ratios, material and equilibrium parameters, and the bifurcation step Ay5 = 0.01 . It can be observed 

from Figure 12(a) that the response of the coupled system has a complete process from chaotic motion 
through periodic motion and chaotic motion to period- 1 motion. At the interval of 0.3 < yff < 0.43, the 
system response enters periodic motion from chaotic motion, then it becomes chaotic motion again, 
and finally it comes into steady-state motion with period- 1 in the range of 0.43 < y? < 0.6. With the 
increase of squeeze film damping //(// = 0.14), the system response changes noticeably and it mainly 
contains the periodic components, such as period- 1, period-3 and period-6 motions, as illustrated in 
Figure 12(c). As the equilibrium parameter a increases, the chaotic components of the system response 
decrease, while the periodic components increase and contain period-2, period-4 and period-8 motions 
with the case of a = 1.6, as shown in Figure 12(d). 

The Poincare maps and phase portraits for different cubic stiffness ratios are displayed in Figure 13. 
When the cubic stiffness ratio is small, i.e., ^ = 0.35, the exhibited motion is period-6 motion, six 
points and circles can be seen in the Poincare map and phase portrait respectively from Figure 13(a). 
As the cubic stiffness ratio becomes larger, the closed circle is decomposed and the points in the 
Poincare map gradually scatter. At y? = 0.42, the system response comes into chaotic motion, and then 
the points of the attractor are decomposed again and finally converge to one point. When y? = 0.48, a 
synchronous motion with period- 1 can be observed. Then the system response comes into period-4 
motion. It is indicated that the components of chaotic motions become wider with the increase of the 
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squeeze film damping. Material properties of the tip and sample and equilibrium coefficient ratio play 
very important role in the nonlinear dynamics of the coupled system. 

Figure 12. Bifurcation diagram of the cubic stiffness ratio y? at different combinations of 

squeeze film damping ratios, material parameters and equilibrium parameters: (a) 
7 = 0.08, 2: = 0.3, a = 1.2; (b) /7 = 0.08, 2 = 0.5, a = 1.2; (c) /7 = 0.14, S = 0.5, 

a = 1.2;(d) 77 = 0.14, 2: = 0.5, 6^ = 1.6. 
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Figure 13. The Poincare maps and phase portraits of different cubic stiffness ratios y?. 
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4.4. Effect of Equilibrium Parameter a 

Equilibrium coefficient ratio a is one of the important parameters for determining the equilibrium 
position of the cantilever tip in AFM and it becomes the key factor to reflect the dynamic responses of 
the tip-sample model. The dynamic behavior of the coupled system depends on the value of a. Figures 
14 and 15 give the bifurcation diagram and largest Lyapunov exponent map of the dynamic system 
with the control parameter of the equilibrium parameter a at = 1. It can be seen from Figure 14 that 
the dynamic responses are very complicated, and the components contain periodic and chaotic motions 
at the interval of 1< a < 2. The corresponding largest Lyapunov exponents are alternately positive and 
negative, as shown in Figure 15. 



Figure 14. Bifurcation diagram of the equilibrium parameter a. 
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To explain the dynamic responses of the system clearly, Figure 16 shows the local bifurcation 
diagram and Poincare maps of the dynamic system at the interval of 1.6 < a < 2. It can be found that 
the system responses exhibit the alternation of periodic and chaotic motions. The system response 
comes into steady-state synchronous motion with period- 1 from chaotic motion, and enters period-2 
motion from period- 1 motion as the equilibrium parameter a increases, and then becomes chaotic 
motion with period-doubling bifurcation. Moreover, at 1.7 < a < 1.9, the system response changes 
between period- 1 and period-2 motions alternately. When a > 1.9, the system response comes into 
steady-state synchronous motion with period- 1. These phenomena indicate that the dynamic responses 
of the coupled system are very complex. The cantilever tip can undergo a period-doubling cascade to 
possible chaos about the original equilibrium. It is demonstrated that, away from the surface, the net 
force on the tip is always in the downward direction and causes the tip to accelerate the sample until it 
passes the key point, where the repulsive force plus the spring force becomes larger than the van der 
Waals force, and then the tip is forced away from the sample. 

Figure 16. Local bifurcation and Poincare maps of equilibrium parameter a . 
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Figure 17 shows the bifurcation diagram of the equilibrium parameter a on the response of TM- 
AFM tip-sample system at the interval of 1.0 < or < 1.8 , the squeeze film damping ratio is taken as 
f] = 0.008 and the bifurcation step is Aa = 0.01 . It can be seen from Figure 17(a) that the response of 
the coupled dynamic system comes into period-2 motion from period- 1 motion and then enters chaotic 
motion with period-doubling bifurcation. With the increase of a , the system responses enters period-2 
motion from chaotic motion, and it subsequently becomes period- 1 motion again when Z = 0.3 and 
r = 2. As illustrated in Figure 17(b), the response of the coupled system has a complete process from 
chaotic motion through periodic motion to chaotic motion with the forms of period-doubling 
bifurcation and anti period-doubling bifurcations at the interval of 1 < 6if < 1.8 when 2 = 0.3 and F = 4. 
With the increases of the material parameter coefficients, as shown in Figure 17(c) and (d), the system 
response changes noticeably and it mainly contains the periodic components, such as period- 1, period- 
2, period-3 and period-6 motions. Meanwhile, the chaotic components of the system response decrease 
and, on the contrary, the periodic components increase. In addition, the chaotic components of the 
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system response shift to the smaller equilibrium parameter. It demonstrates that the components of 
chaotic motions become wider as the amplitude of external forcing term increases, and the material 
properties of the tip and sample and equilibrium coefficient ratio affect the nonlinear dynamics of the 
coupled system. 

Figure 17. Bifurcation diagram of the equilibrium parameter a at different combined 
material parameters and external forcing terms: (a) E = 0.3, r = 2;(b) Z = 0.3, r = 4; 

(c) z = o.5, r = 2;(d) s^o.s, r = 4. 




To illustrate the various motions. Figure 18 shows the nonlinear characteristics of the coupled 
system with the plots of the Poincare maps and phase portraits for different equilibrium parameters a 
at different conditions. The motion of the coupled system changes between periodic and chaotic 
motions alternately. At a = 1.2, the motion with period- 1 represented by a point in the Poincare maps 
and characterized by a close curve in phase portraits is shown in Figure 18(a). As illustrated in Figure 
18(b), the system response comes into period-6 motion at a = 1.42 from synchronous motion with 
period- 1 at a = 1.2, as displayed in Figure 18 (a), then leaves period-6 motion and enters chaotic 
motion at a = 1.6, which can be seen from Figure 18(c). The strange attractor has a fractal structure 
and the corresponding largest Lyapunov exponent is positive. With the increase of the equilibrium 
parameter coefficient, as shown in Figure 18(d), the system response becomes periodic motion from 
chaotic motion again, and one can find the period-3 motion marked by three isolated points in Poincare 
map and three circles in phase portrait at a = \.lA. It is indicated that the components of chaotic 
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motions of the coupled system increases obviously with the increase of the amplitude of the force 
term. In general, the effect of equilibrium parameter on the system response should be considered for 
the design of the TM-AFM. 

Figure 18. The Poincare maps and phase portraits of different equilibrium parameters a . 

(a) a^\2 (b) 6^ = 1.42 (c) a = 1.6 (d) = 1.74 




5. Conclusions 



The cantilever in tapping mode Atomic force microscope (TM-AFM) is one of crucial components 
and so it is very important to carry out a thorough dynamic analysis for the cantilever to further 
enhance the performance of the TM-AFM. In this paper, the cantilever- sample interaction and a 
harmonically forcing term in the TM-AFM have been considered. Numerical simulations have been 
used to investigate the nonlinear behaviors between the tip and the sample. The chaotic behavior 
appears to be generated via various system parameters including the amplitude of the external forcing 
term, equilibrium parameter, squeeze film damping and material property. The dynamic system 
responses display very rich nonlinear dynamic characteristics under the effects of these parameters and 
show an alternate changing process among periodic motion, quasi-periodic motion and chaotic motion. 
The component of the chaotic motion in the system response increases with the increase of the 
amplitude of the external excitation, but will be weaken by the squeeze film damping. In addition, the 
increase of the equilibrium parameter coefficient would strengthen the component of the chaotic 
motion but weaken the component of the periodic motion increases in the system response. It is 
indicated that the external excitation and squeeze film damping are important for the design of 
dynamic TM-AFM. It is of significance to understand the cantilever dynamics in air/liquid and 
promote the development of the next generation of AFMs. 
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